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Abstract: We present spectra of highly excited D and D s mesons up to around 3.8 GeV 
determined using dynamical lattice QCD. We employ novel computational techniques and 
the variational method with a large basis of carefully constructed operators in order to 
extract and reliably identify the continuum spin of an extensive set of excited states. 
These include states with high spin and states identified as having an explicit gluonic 
contribution. Calculations were performed on two volumes, both with a pion mass of 
approximately 400 MeV, achieving a high statistical precision for both ground and excited 
states. We discuss our results in light of experimental observations, comment on the 
phenomenological implications and identify the lightest 'supermultiplet' of hybrid mesons 
in each sector. 
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1 Introduction 

The spectra of open-charm mesons contain a number of experimentally well-established 
states [1]: eight "charm- light" D mesons with isospin 1=1/2 and six "charm-strange" 
D s mesons with non-zero strangeness and zero isospin. Many of the states whose spin, J, 
and parity, P, are well determined fit the pattern expected from model predictions but a 
subset of these, as well as additional states which require experimental confirmation, do not 
and remain unexplained. Suggestions such as hybrid mesons (in which the gluonic field is 
excited), molecular mesons and tetraquarks (with four valence quarks) have been proposed, 
but no clear picture has emerged. In particular, the masses and widths of the enigmatic 
D*g(2317) =1: and Z?< j i(2460) ± are lighter and narrower than in quark models [2, 3]. Since 
the discovery of these resonances at BABAR [4] and CLEO [5] respectively, there has been 
much debate about their nature and the situation is yet to be resolved. 

A simple picture of a meson describes it as a bound state of a quark and an antiquark. 
If the meson is an eigenstate of charge-conjugation with quantum number C, some J PC 
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combinations are inaccessible in this description; such states are termed 'exotic'. The 
charmonium system provides an arena for experiments such as BESIII and PANDA to 
search for these exotics, with the observation of such a state being a smoking gun for 
physics beyond quark models. A recent study of charmonium by the Hadron Spectrum 
Collaboration found states with exotic quantum numbers and showed that these, along with 
other non-exotic states, could be identified as hybrid mesons [6]. On the other hand, open- 
charm mesons are not eigenstates of charge conjugation and so there are no such exotic 
quantum number combinations. This increases the difficulty of experimentally identifying 
a multiquark or hybrid state in the open-charm sector. However, hybrid states are still 
expected to be present in the spectrum and the system provides an important probe of the 
dynamics of excited glue in another environment. 

A full understanding of QCD and its spectrum requires ab initio calculations to com- 
pare against and inform experimental searches, and to discriminate between models to 
glean the relevant physics. Spectroscopic calculations on the lattice have a long history 
and an accurate determination of the low-lying experimental spectrum is considered an 
acid test of the precision and accuracy of the method. Current lattice QCD calculations 
have reached this benchmark as shown by the extremely precise calculations of the so-called 
"gold plated" hadron spectrum [7]. 

The extraction of higher-lying spectroscopic states on the lattice is more problematic. 
In lattice calculations space-time is discretised on a hypercubic grid of points and its 
volume reduced to that of a few fermi. The breaking of Lorentz invariance reduces the full 
continuum rotation symmetry to that of the octahedral group, Oh', this means that states at 
rest are no longer classified according to their spin, J, but by the irreducible representations 
(irreps) of the octahederal group. Since there are only five irreps 1 of Oh for each parity 
and an infinite number of spins, each irrep contains an infinite number of spins and the 
various components of states with J > 2 are distributed across several irreps; this makes the 
task of disentangling the spin of states difficult. The Hadron Spectrum Collaboration has 
addressed the problem of this reduced symmetry through the use of large bases of carefully 
constructed interpolating operators. In combination with dynamical anisotropic lattices [8, 
9] and the distillation technique [10] , this methodology has recently been used successfully 
to extract highly excited spectra and to reliably identify the continuum quantum numbers 
of high spin states for both light and heavy mesons [6, 11-13], as well as for baryons [14- 
17]. We will use the same methodology in this study and summarise the important points 
below. 

Previous calculations of the open-charm meson spectra have mainly focused on deter- 
mining the lowest-lying S'-wave [(0, 1)~] and P-wave [(0, 1, 2) + ] states (see Refs. [18-20] for 
some recent investigations). Only more recently have calculations begun to explore higher 
up in the spectrum [21-24] and we discuss these studies in Section 5.2. In this paper we 
extract D and D s meson excitation spectra at a single lattice spacing with a pion mass of 
approximately 400 MeV and we reliably identify the continuum quantum numbers of the 
extracted states. The results presented here are the most extensive determination of the 

1 here we are only considering the single-cover representations relevant for mesons 
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open-charm spectra from lattice QCD to date, include all J p combinations with J < 4 
and consider states with a hybrid nature. Preliminary results from this work have been 
presented in conference proceedings [25]. 

In addition, this work lays a foundation for more detailed investigations of resonances 
above strong-decay threshold in both charmonium and the open-charm sector, including 
the enigmatic new states. These studies will consider scattering involving D and D s mesons 
for which the determination of masses and dispersion relations discussed in this paper is an 
important prerequisite. We can compute the multi-meson energy levels in a finite volume 
by supplementing our operators with those designed to efficiently create such states [26, 27]. 
The resulting denser spectrum can be analysed, at least in the case of elastic scattering, us- 
ing the Liischer methodology [28] and its extensions, and the mass and width of resonances 
determined. 

2 Lattice discretisation and interpolating operators 

The calculations presented here make use of ensembles generated by the Hadron Spectrum 
Collaboration; for more details see Refs. [8, 9]. These gauge field configurations include the 
dynamics of two mass-degenerate light quarks and a strange quark, while the charm-quark 
field is quenched. We use an anisotropic discretisation in which the spatial lattice spacing, 
a s , and the temporal lattice spacing, at, are distinct and related via £ = a s /at ~ 3.5. 
This ensures that we have a^m c <C 1, where m c is the bare charm-quark mass, and the 
standard relativistic formulation of fermions can be used to study states in which the charm 
quark four-momentum is closely aligned with the temporal axis. It is worth noting that in 
these computations a s m c is also less than unity. The anisotropy, £, is discussed further in 
Section 4. 

2.1 Gauge and light-quark fields 

In the gauge sector, a Symanzik-improved anisotropic action with tree-level tadpole-improved 
coefficients is used. For all the quark fields, a tree-level tadpole-improved Sheikholeslami- 
Wohlert (clover) anisotropic action with stout-smeared spatial links [29, 30] is employed. In 
Refs. [8, 9] the parameters of these actions were determined in order to study light-hadron 
spectroscopy using dynamical light and strange quarks. 

In lattice computations, all quantities are calculated in terms of the lattice spacing 
and, to make contact with experiment and present results in physical units, a scale-setting 
is introduced. Here we follow Ref. [9], using the ratio of the O-baryon mass measured on 
these ensembles, atm^i = 0.2951(22) [15], to the experimental mass, mn = 1672.45(29) 
MeV [1], to determine a t . This yields a^ 1 = 5.67(4) GeV and a spatial lattice spacing 
about 3.5 times larger, a s ~ 0.12 fm. We perform calculations on two lattice volumes 
corresponding to spatial extents of about 1.9 fm and 2.9 fm. Table 1 summarises the two 
ensembles used in this study, with full details given in Refs. [8, 9]. 
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Lattice size m^/MeV -/V c f gs iVtsrcs 



16 3 x 128 391 96 128 64 

24 3 x 128 391 553 16 162 



Table 1. The gauge-field ensembles and quark propagators used in this work on two lattice volumes, 
(L/a s ) 3 x (T/at), where L and T are respectively the spatial and temporal extents of the lattice. 
The number of gauge-field configurations used, N c [ gs , and the number of time-sources for quark 
propagators per configuration, N tS rcs, are shown; N vecs refers to the number of eigenvectors used in 
the distillation method [10]. 

2.2 The charm-quark action 

Lattice QCD has a long-standing problem when trying to perform calculations involving 
heavy quarks moving relativistically. One requires very fine lattices because the usual 
lattice actions lead to large discretisation artefacts when the lattice spacing a > m^ 1 , where 
mo is the bare quark mass. In theory there is a simple solution: perform calculations on 
a much finer lattice. But in practice, this is computationally expensive if a large enough 
spatial volume is maintained. Recently, anisotropic lattices have been used to circumvent 
the problem by employing a temporal lattice spacing smaller than that in the spatial 
direction. This enables simulations with at mo <C 1, allowing correlation functions to be 
studied in detail and discretisation artefacts of heavy quarks to be kept under control, while 
maintaining the spatial directions coarse enough to keep the computational cost reasonable. 

In Ref. [6] , a valence charm quark was introduced and used to study the charmonium 
spectrum. The bare mass of the charm quark was non-perturbatively determined by re- 
quiring the ratio of the r] c mass to the O-baryon mass to take its experimental value. The 
same action was used for the charm-quark field as for the light and strange quarks, except 
the parameter giving relative weights to spatial and temporal finite differences was tuned to 
give a relativistic dispersion relation, consistent with £ = a s /at = 3.5, for the r/ c meson at 
low momentum. In this study we use the same actions and set of parameters. In Section 4, 
we investigate D and D s mesons with non-zero momentum and find that they satisfy a 
relativistic dispersion relation, again consistent with £ = 3.5. 

2.3 Operator construction 

Spectral information can be obtained from Euclidean two-point correlation functions de- 
termined by lattice QCD computations, 

C ij (t) = (Q\O i (t)O}(0)\0), (2.1) 

where 0^(0) anc i 0(t) are respectively the source and sink interpolating fields (creation 
and annihilation operators). By inserting a complete set of eigenstates of the Hamiltonian, 
the correlation function becomes a sum over all states with the same quantum numbers of 
the interpolating fields, 
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Table 2. Gamma matrix combinations and naming scheme. 
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Table 3. Distribution of continuum spins across the irreps of Oh, A(dim), where dim is the 
dimension of the irrep. 



with the discrete nature of the spectrum arising from the finite volume of the lattice. 
The vacuum-state matrix elements Zf = (n|O||0), or overlaps, are used to identify the 
continuum spin of the extracted states as discussed in Section 3.1 and to probe their 
structure. 

In order to maximise the spectral information we can extract from correlation functions 
we employ the following strategies. Firstly, we construct a large basis of operators with a 
variety of spin and spatial structures that allows us to explore all J p combinations with 
J < 4. In general, these operators take the form of a quark bilinear with an operator 
insertion T, 

O ff =tpfT V?/' = i>f T^i^Bj ■ ■ ■ ipf , (2.3) 

where /,/' label the quark flavours, T is a gamma matrix combination taken from Table 2 
and the projection onto definite momentum and spatial indices are suppressed for clarity. 
The 'forward-backward' covariant derivative is defined as D = As described 

in Refs. [11, 12], the derivatives and vector-like gamma matrices are first expressed in a 
circular basis so that they transform as J = 1. Using the usual SO(3) Clebsch-Gordan 
coefficients they can then be combined to form operators, O j , with any desired continuum 
J p . In order to use these operators in lattice calculations they are subduced into the 
relevant irrep(s) of Oh [12] and Table 3 shows how the continuum spins are distributed 
across these irreps. Operators subduced into irrep A and irrep row A from spin J are 
denoted by C>j^ A and our operator naming scheme is explained in Ref. [12]. The number 
of subduced operators employed in each irrep is shown in Table 4. 

The operator in the bilinear, T, is constructed to have well-defined transformation 
properties under transposition, T T = ±I\ We denote the eigenvalue of Y under this 
transformation as Cr = ±1. If the quark and antiquark fields have the same flavour, 
i.e. f = /', this corresponds to the charge-conjugation quantum number, C. Similarly, if 
/ 7^ /' but / and f are degenerate in mass, this is related to a quantum number which 
is a generalisation of C, e.g. G-parity relevant for mesons consisting of degenerate up and 
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Table 4. The number of operators used in each lattice irrep, A p . All combinations of gamma 
matrices and up to three derivatives are included. 

down quarks. For D and D s mesons, / = c, /' G {u, d, s} and there is no mass degeneracy 
between the quark and antiquark fields in our operators. As a result, states created by 
operators with opposite Cr can mix to form an eigenstate of the QCD Hamiltonian. This 
effect will be considered in more detail in Section 3. 

Secondly, we employ a now well-established quark smearing procedure known as dis- 
tillation [10]. This procedure consists of applying a smoothing function to the quark fields, 
reducing the contamination from noisy UV modes that do not make a significant contri- 
bution to the low-energy physics we wish to extract. The distillation operator is defined 
via, 

^(t)=J2f^ n) )4 n \t)vf\t), (2.4) 

n=l 

where (t) are the eigenvectors, sorted by eigenvalue \( n \ of the three-dimensional gauge- 
covariant lattice Laplacian evaluated on the background of spatial gauge fields of timeslice 
t. f(\( n ') is a smearing profile and in this study it is set to unity, although other smearing 
profiles could prove to be beneficial. Apart from the removal of short-distance modes, 
distillation has other significant benefits. One advantage arises from the outer product 
structure of the distillation operator, □; the correlation function can be factorised into per- 
ambulators and elementals with the former containing information on the propagation of 
the quarks and the latter encapsulating the momentum and structure of the state created. 
This factorisation allows the perambulators to be stored and so reduces the computational 
cost of the calculation. Once the linear systems arising from the Dirac equation on a given 
gauge background have been solved for as many time-sources as required, correlation func- 
tions involving any elementals can be computed without further inversions. This enables 
us to perform calculations with large bases of operators having a variety of spin and spa- 
tial structures. In addition, distillation enables the efficient computation of annihilation 
contributions and of correlation functions involving multi-hadron operators. 

Using distillation, for each irrep we compute two-point correlation functions for all 
possible combinations of our operators. Our methodology for analysing these correlators 
is described in the following section. 
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3 Correlator analysis 



For each lattice irrep and flavour sector (D and D s ) we compute an N x N matrix of 
correlation functions, 

C ij (t) = (0\O i {t)o](p)\0), (3.1) 

where i,j label the operators in the basis and N is the number of operators in the irrep, 
given in Table 4. We then apply the variational method [28, 31], now commonly used in 
lattice QCD calculations of spectra, which finds the optimal extraction, in a variational 
sense, of energies in a given channel. The details of our particular implementation along 
with various systematic tests are given in Ref. [12] and tests in the charmonium sector have 
been presented in Ref. [6]; here we will briefly summarise the important points. Practically, 
the procedure is realised via the solution of a generalised eigenvalue problem, 

Cyflv] = \ n (t, t )C l3 (t )v], (3.2) 

where the reference time-slice to must be chosen appropriately as described in Refs. [12, 
32]. Solving this for each t yields eigenvalues, A n (t, to), and eigenvectors, v™, where n = 
1,2, . . . , N . The method is very powerful for extracting excited-state energies; even nearly 
degenerate states can be distinguished because such states in general have different patterns 
of operator-state overlaps corresponding to distinct orthogonal eigenvectors. 

The eigenvalues, known as principal correlators, allow access to the energies, E n , via 
their dependence on (t — to); for sufficiently large times they are proportional to e ~ En ^~ to \ 
The eigenvectors are related to the overlaps, Zf, which give information on the structure of 
a state and, as will be described below, enable us to identify the continuum spin of a state. 
For insufficiently large times, there will be some contamination from excited states in the 
principal correlators. In order to 'mop up' this contamination and enable us to stabilise 
the fits by considering a larger temporal extent, we use a fitting function containing two 
exponentials of the form, 

A n (t) = (1 - A n )e~ mn(t - to) + A n e- m " ( '-' o) , (3.3) 

where the fit parameters are m n , m' n and A n . Empirically we find that the contribution 
from the second exponential decreases rapidly as we increase to- Our best estimate of the 
mass of the state is m n and we do not further consider the quantities m' n and A n . 

The combination of distillation and variational analysis of correlator matrices with 
large bases of subduced operators has proved very effective at extracting spectra of excited 
and high-spin light mesons [11-13], charmonia [6] and baryons [15, 16]. 

3.1 Spin identification 

In principle, to determine the continuum spin, J, of states extracted in lattice QCD cal- 
culations, the spectrum can be computed at finer and finer lattice spacings and then ex- 
trapolated to the continuum limit. Exact degeneracies will then emerge between states 
of the same spin subduced into different irreps of Oh', these can be matched up between 
irreps enabling the spin to be determined from the pattern of degeneracies. However, there 



-7- 



are a number of problems with this approach. The high computational cost of performing 
calculations at multiple fine lattice spacings means that this is not currently feasible for 
computations of excited states. More problematically, there can be physical degeneracies 
or near-degeneracies in the dense high-lying spectrum; the question then arises as to how 
to distinguish these, with a finite statistical precision, from artefacts due to a finite lattice 
spacing. It would be useful to have a spin identification scheme that can be used to reli- 
ably determine the spin at a single lattice spacing. As we will see in Section 5, the dense 
spectrum of excited states would be impossible to disentangle without more information 
beyond solely the masses of states. This is particularly relevant in the D and D s meson 
sectors where a lack of charge-conjugation symmetry makes the spectrum more dense than, 
for example, the charmonium sector. 

The Hadron Spectrum Collaboration has recently developed a scheme [11, 12] to iden- 
tify the continuum spin of an extracted state at a single lattice spacing by using information 
from the operator-state overlaps defined in Section 2.3. Further details and a demonstration 
of the efficacy of this approach are presented in Ref. [12] and here we will summarise the 
method; an application to the charmonium sector is given in Ref. [6]. The subduced opera- 
tors, respect the symmetry of the lattice but we find that they also carry a 'memory' 
of the continuum spin, J, from which they were subduced. If full rotational symmetry 
reasonably describes transformations of the smeared fields on hadronic scales, we would 
expect that an operator subduced from a continuum spin J will overlap predominantly 
onto states of continuum spin J [33]. 

Figs. 1 and 2 show that these statements are justified at the level of the correlation 
matrix and the overlaps respectively. In Fig. 1 we show the normalised correlation matrix 
at time-slice 5 for the irrep in the charm-light (D) sector. The matrix is approximately 
block diagonal in J, suggesting that, for example, an operator subduced from J = 1 has 
little contribution from the J = 3 and J = 4 states. As discussed in Section 2.3, charm-light 
and charm-strange mesons are not eigenstates of charge conjugation. However, our operator 
insertions have a definite symmetry under transposition (equivalent to charge conjugation 
of the operator followed by the interchange of the quark and antiquark flavours) and we 
use this to label the operators, Cr = ±. This would correspond to a generalisation of 
charge-conjugation parity (analogous to G-parity) if the charm and light/strange quarks 
were degenerate, i.e. if we had SU (4) symmetry. The lack of such a symmetry is shown in 
Fig. 1 where significant overlap between Cr = + and Cr = — operators is observed. 

In Fig. 2 we show the overlaps of a selection of operators onto a selection of low-lying 
states in the irrep of the charm-strange (D s ) sector. Every state clearly shows a strong 
preference to overlap onto only operators subduced from a particular spin. The lack of a 
generalised charge-conjugation symmetry is shown by the observation that the J = 1 state 
has strong overlap onto both Cr = + and Cr = — operators. The subset of overlaps we 
show here is representative of the full operator basis - the observation that each state only 
overlaps strongly onto operators subduced from a particular spin is repeated across our 
entire set of spectra. 

For J > 2, we can also match up the states across different irreps of Oh by comparing 
the absolute values of the overlaps. Up to discretisation effects, these values should be 
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Figure 1. The normalised correlation matrix, Cy/ y/CuCjj, on time-slice five for the T± irrep in 
the charm-light (D) sector on the 24 3 volume. The operators are ordered according to continuum 
spin, J, as indicated by the labelling. For each value of J the operators are then ordered according 
to their symmetry, Cr = ±, as described in the text. The correlation matrix is observed to be 
approximately block diagonal in spin. 



consistent for states subduced from the same continuum state with spin J. In Fig. 3 we plot 
a selection of overlaps for states conjectured to have J p = 2~,3~ and 4~, demonstrating 
there is good agreement between the overlaps extracted in different irreps. This observation 
is repeated across our entire spectra though, as we expect, there are slight deviations from 
equality due to discretisation effects. To quote masses for J > 2 states we perform joint 
fits to the principal correlators in each irrep as described in Ref. [12]. 

In summary, we have demonstrated that the method for identifying the continuum 
spin of extracted states that was developed in Refs. [11, 12] is applicable to studies of 
charm-light and charm-strange mesons. We will use this approach to identify the J of the 
extracted states that we present in Section 5. 

4 Dispersion relations 

One purpose of this investigation is to ensure we can study heavy-light mesons reliably 
on these anisotropic lattices to enable subsequent calculations involving scattering of these 
states. This is essential if the nature of the charmonium and charmed resonances above the 
open-charm threshold is to be understood. To determine scattering properties of mesons 
using lattice QCD requires the computation of the full spectrum of energy eigenstates, 
including those made up predominantly of more than one meson. For these determinations 
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Figure 2. The overlaps, Z, of a selection of operators onto a selection of states in the charm- 
strange (D s ) T x + irrep on the 24 3 volume. In this plot the overlaps have been normalised so that 
the largest value for an operator across all states is equal to unity. The error bars indicate the 
one sigma statistical uncertainty. For each state, the operators are coloured from left to right 
as: pink (a x D [ f =l f J=l \ red (h x D [ f =0 )^= l \ blue (b x x D [ f =2 f J = 3 \ red (p x D^)^, pink 
( P xDf i3=2 J=2 )^\ light blue (/)xD^ 2 J=2 )[«] and gold (pxD^y^)^]. The first three 
operator insertions have negative symmetry (Cr = — ) and the last four have positive symmetry 
(Cr = +) as explained in the text. States are labelled by their mass in temporal lattice units and 
the continuum spin of the dominant operators. 'Hyb' refers to a state which has relatively strong 
overlap with operators that are proportional to the field strength tensor, the commutator of two 
gauge-covariant derivatives. 
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Figure 3. A selection of absolute ^-values for the lightest states conjectured to have J = 2 (left 
pane), J = 3 (middle pane), J = 4 (right pane) in charm-strange (D s ) irreps, A~, on the 24 3 
volume. For each J we show the values of Z for the same operator subduced into different irreps. 
The vertical size of a box corresponds to the one-sigma statistical uncertainty on either side of the 
mean. 
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Figure 4. Squared energies as a function of |n| 2 for the lightest S-wave and P-wave states in both 
the D (left panel) and D s (right panel) sectors. The points correspond to calculated energies and 
the statistical uncertainties are smaller than the size of the points. The lines are dispersion relations 
of the form of Eq. 4.1, as described in the text, with m fixed to the rest energy of the appropriate 
meson and £ fixed to 3.46 for the charm-light mesons and 3.45 for the charm-strange mesons. The 
measured energies are seen to be well described by the dispersion relations. 

to be reliable, the dispersion relations for the mesons produced by decays of the resonances 
should be relativistic up to the appropriate momenta set by the decay kinematics. Heavy 
quarks bring extra complications to these calculations but the anisotropic discretisation 
we use here may reduce the impact of the large mass scale, mg. For mesons moving with 
modest momenta, the heavy quark four-momentum will predominantly be aligned with the 
temporal direction, the most finely discretised lattice direction. 

One well-studied artefact of heavy quarks on the lattice is that when awiQ is 0(1) 
the dispersion relation for mesons appears non-relativistic and the kinetic and rest masses 
of the mesons differ. Heavy-light and heavy- heavy states are effected differently. On the 
anisotropic lattice, the relativistic dispersion relation is most naturally written, 

{a t E? = (a t m) 2 + QV \a s \p\) 2 , (4.1) 

where p = ^-n follows from the periodic boundary conditions of the spatial directions. 
Here n = (n x ,n y ,n z ) with rii taking integer values. To determine whether the anisotropic 
lattice has controlled the size of these artefacts, we determined the dispersion relations for 
a number of mesons in the D and D s spectra. Recall that the parameters in the action 
for the charm quark are determined such that the r/ c meson has a relativistic dispersion 
relation consistent with £ = 3.5; we obtained £„ c = 3.50(2) [6]. This compares with an 
anisotropy measured from the pion dispersion relation of £ n = 3.444(6) [26]. If there was 
a problem with large discretisation artefacts, one symptom could be that the £ measured 
from the D or D s dispersion relations would be different from the £ determined from the 
pion or the r/ c . 

To calculate dispersion relations for a subset of D and D s mesons, we use two sets of 
operators. For mesons at rest we use the same basis of operators as described in Section 2.3, 
but for mesons of non-zero momenta we follow Ref. [34] and use operators constructed in 
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a helicity basis (with up to two derivatives) and then subduced into the relevant irreps. In 
Fig. 4 we show the squared energies as a function of \n\ 2 for the low- lying states in both the 
charm-light and charm-strange sectors. The points correspond to calculated energies and, 
as usual, if the state subduces into more than one irrep, these energies are from joint fits to 
the relevant principal correlators. We determine the anisotropics in the D and D s sectors, 
£d and £d s , by fitting Eq. 4.1 to the ground state (0~) D and D s energies with \n\ 2 < 4. 
This gives good fits with xV^dof ~ 1 yielding £d = 3.46(1) and £d s = 3.45(1). The lines 
in Fig. 4 are of the form of Eq. 4.1 with m fixed to the rest energy of the relevant meson 
and £ fixed to our measured 0~ D or D s anisotropy as appropriate. We do not see much 
deviation of the higher energy states from the dispersion relations using these anisotropics. 
This suggests that we are correctly describing relativistic states up to at least \n\ 2 = 4 
with an anisotropy not differing greatly from the measured pion and rj c anisotropies. This 
is crucial for future calculations of scattering involving open-charm mesons. 

5 Results 

In this section we present the results of our study. We first compare the energy levels 
extracted in each lattice irrep across the two volumes, then give final spin-identified spectra 
and comment on other lattice investigations. We discuss the implications of our results 
and compare to experiment in Section 6. 

5.1 Results by lattice irrep and volume comparison 

The results of the variational analysis for each lattice irrep in the charm-light (D) and 
charm-strange (D s ) sectors are shown in Figs. 5 and 6 respectively. The vertical size of 
the boxes indicates the one sigma statistical uncertainty on either side of the mean while 
ellipses indicate that additional states may be present in the variational analysis in this 
energy region but these were not robust. The colour of a box represents the continuum spin 
assignment from the spin identification scheme discussed in Section 3.1: states identified 
as J = are coloured black, J = 1 are red, J = 2 are green, J = 3 are blue and J = 4 are 
orange. 

Components of J > 2 states are distributed across more than one lattice irrep and 
levels appear at similar masses, up to discretisation effects, in different irreps. The vari- 
ous components of the same continuum state are matched by comparing the overlaps as 
described in Section 3.1. From Figs. 5 and 6 it is clear that there are no significant dif- 
ferences between energies of the same continuum state subduced into different irreps. The 
extremely dense spectra of states appearing at and above a^m ~ 0.52 in the figures would 
be impossible to disentangle using only the extracted masses, emphasising the value of our 
spin identification scheme. 

In general, throughout our spectra we find no significant volume dependence. However, 
there are a couple of exceptions: for the lightest + and 1 + , determined with high statistical 
precision, we find a 2a discrepancy between the two volumes. These states lie precariously 
close to thresholds and therefore mixing with multi-meson states may be important and 
could be the cause of the observed volume dependence. Our operator bases do not include 
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Figure 5. The extracted spectrum of states in the D meson sector labelled by irrep A p . For each 
irrep results from both the 24 3 and 16 3 volumes are shown side by side. The vertical size of each box 
gives the one sigma statistical uncertainty on either side of the mean and the box colour refers to 
the continuum spin assignment as described in the text. The light cyan boxes represent states that 
were not very well determined in the variational analysis; ellipses indicate that additional states 
may be present in these energy regions but were not robustly determined. 
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Figure 6. As Fig. 5 but for the D s meson sector. 
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any operators that 'look like' two-mesons and so, as discussed in Refs. [12, 27], we do not 
expect to be able to reliably extract multi- meson energy levels; a conservative approach is 
to suggest that our mass values are accurate only up to the hadronic width. Note that in the 
charm-strange sector these states correspond to the enigmatic D* (2317) =l: and -D s i(2460) =t 
which have been suggested to be molecular states of two mesons or have their properties 
modified by interaction with the nearby thresholds. To investigate this in more detail, in 
future work we will include multi-meson operators in our bases. This will result in a denser 
spectrum that can be analysed along the lines described in Refs. [26, 27]. Ref. [22] has 
recently performed some investigations of the lightest + and 1 + charm- light resonances in 
simulations with dynamical up and down quarks (strange-quark degrees of freedom were 
not included). 

With these possible exceptions we find no significant volume dependence and the pat- 
tern of states across irreps is consistent with components subduced from states of definite 
continuum spin (multi-hadron states would show a different pattern). We therefore gener- 
ally see no clear evidence for multi-hadron states in our spectra. 

5.2 Final spin-identified spectra 

Our final results are the spin-identified spectra from the 24 3 volume - these were obtained 
with higher statistics than those from the 16 3 volume and we can reliably extract and 
identify the spin of a larger number of states. In addition, the large volume means that 
any possible finite- volume effects will be less important. We show our spin-identified spectra 
for charm-light (D) and charm-strange (D s ) mesons in Figs. 7 and 8 respectively; here we 
only show well determined states for which the spin could be reliably determined. We 
show the calculated (experimental) masses with half the calculated (experimental) r] c mass 
subtracted in order to reduce the systematic error from any imprecision in the tuning of 
the bare charm-quark mass (see Section 2). In Fig. 7, the dashed lines correspond to 
the lowest non-interacting Dir and D S K thresholds using our calculated masses (coarse 
green dashing) and using experimental masses (fine black dashing). The dashed lines in 
Fig. 8 correspond to the lowest non-interacting DK threshold using our calculated masses 
(coarse green dashing) and using experimental masses (fine black dashing). These results 
are tabulated in Appendix A. 

In the next section we give some interpretation of our results and compare with the 
experimental situation but first we discuss other lattice computations of excited charm-light 
and charm-strange spectra. Most lattice studies have concentrated on the lowest-lying S 
and P-wave states below threshold (see Refs. [18-20] for some recent studies) and only 
in the last couple of years have calculations begun to explore states higher up in the 
spectrum [21-24]. These computations, like ours, have various systematic uncertainties 
which have not been fully accounted for and so a direct quantitative comparison is not 
appropriate. Depending on the study, these systematics come from having a finite lattice 
spacing, working in a finite volume, having light quarks with unphysical masses and not 
taking into account the finite width of states. 

In Ref. [21] the focus is on the lowest-lying S and P-wave multiplets though the first 
excited S-wave multiplet of D s mesons is also presented. Dynamical (Nf = 2+1) ensembles 
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Figure 7. The D meson spectrum up to around 3.8 GeV labelled by J p . The green boxes are 
our calculated masses on these ensembles with M n ~ 400 MeV, while the black boxes correspond 
to experimental masses of neutral D mesons from the PDG summary tables [1]. We present the 
calculated (experimental) masses with half the calculated (experimental) rj c mass subtracted to 
reduce the uncertainty from tuning the bare charm-quark mass. The vertical size of each box 
indicates the one sigma statistical uncertainty on either side of the mean. The dashed lines show the 
lowest non-interacting Dn and D S K thresholds using our measured masses (coarse green dashing) 
and experimental masses (fine black dashing). 
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Figure 8. As Fig. 7 but for the D s meson sector. The dashed lines indicate the lowest non- 
interacting DK threshold using our measured masses (coarse green dashing) and using experimental 
masses (fine black dashing). 
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were used on a single volume for a number of different pion masses from 702 MeV down to 
156 MeV; the Fermilab method [35] was used for the charm quark. Ref. [22] presents results 
from an Nf = 2 calculation (i.e. with no strange-quark fields) of the charm-light spectrum 
on a single volume with M n = 266 MeV using the Fermilab method for the charm quark 
with the correlation functions computed using the distillation technique. The lowest-lying 
± , 1 ± , 2 ± states are obtained along with some excited and 1 states, and the lightest 
resonances in the + and 1 + channels are investigated. 

Refs. [23, 24] present preliminary results from dynamical computations of charm meson 
spectra for J < 3 at the SU(3) point (strange quarks degenerate with the up and down 
quarks) with M n = 442 MeV using the SLiNC action on a single volume. In all J p channels 
they obtain a ground state and a first excited state while for CP and + they also calculate 
a second excited state. The treatment of the strange quarks in these calculations is very 
different from our approach. The preliminary computations have SU (3) flavour symmetry 
with Mk = M n = 442 MeV and therefore it is not straightforward to directly compare 
with our results. 

The operator bases used in these studies are significantly smaller than ours, they do not 
include as wide a range of spatial structures and in general not all the relevant lattice irreps 
are considered. This means that a more limited number of states can be reliably extracted 
and makes a robust identification of the spin of excited and high-spin states difficult. In 
addition, as discussed in the following section, our inclusion of operators proportional to 
the commutator of gauge-covariant derivatives allows us to explore the hybrid content of 
our spectra. 

6 Interpretation of results and hybrid phenomenology 

In this section we give some interpretation of our spectra. We then discuss and compare 
our results to the experimental situation and explore the mixing between spin-singlet and 
spin-triplet axial-vector, tensor and vector hybrid states. 

Our final results are shown in Figs. 7 and 8 for the charm-light and charm-strange 
mesons respectively, with both spectra showing a similar pattern of states. Most of the 
states fit into the n 2S+1 Lj classification expected by quark potential models, where S is 
the spin of the quark-antiquark pair, L is the relative orbital angular momentum, n is the 
radial quantum number and J is the total spin of the meson. Following Refs. [36, 37], the 
operator-state overlaps, Z, are used to identify the structure of extracted states and make 
2S+1 Lj assignments. 

In both spectra, in the negative-parity sector we find a ground state S"-wave pair 
[CP, 1 _ ] along with a first excitation ~ 700 MeV higher. We see a full D-wave [(1, 2, 3) _ , 2~] 
set at A^ ~ 1400 MeV, where Ae is the mass with half the r) c mass subtracted. We observe 
the second excitation of the S'-wave at A^; ~ 1900 MeV and what appear to be parts of 
an excited D-wave set and parts of a G-wave [(3, 4, 5)~, 4~] at A^; ~ 2000 MeV. The 
5~ state needed to complete the G-wave multiplet is not obtained in our calculation; to 
reliably extract and identify a spin-five state we would need to include additional operators 
that overlap with such a state in the continuum, requiring operators with at least four 
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Figure 9. The negative-parity D (left panel) and D s (right panel) meson spectra showing only 
channels where we identify hybrid candidates. The red boxes are identified as states belonging to 
the lightest hybrid supermultiplet as discussed in the text and other notation is as in Figs. 7 and 8. 



derivatives. We discuss below the supernumerary states appearing at energies between the 
D-wave and second excitation of the S-wave. 

Moving to the positive-parity sector, we observe a full P-wave set [(0, 1, 2) + , 1 + ] around 
the D S K threshold in the charm-light spectrum and around the DK threshold in the charm- 
strange spectrum. We find a first excited P-wave set ~ 600 MeV higher in energy. A full 
F-wave set [(2, 3, 4) + , 3 + ] is present at Ae ~ 1700 MeV in both spectra. We discuss below 
the positive-parity states appearing around = 1900 MeV. 

In both the charm-light and the charm-strange spectra, in the negative parity sector 
we observe four states ~ 1200 MeV above the lowest S*-wave states that do not appear 
to fit into the pattern expected by quark models; they are highlighted in red in Fig. 9. 
Unlike other states, these have relatively strong overlap with operators proportional to the 
field strength tensor on smeared gauge fields, suggesting that the hadronic-scale gluonic 
field plays an important role, and so following Ref. [12, 37] we identify these as the lightest 
"supermultiplet" of hybrid mesons. As explained in Ref. [37], states within a given su- 
permultiplet should have similar operator-state overlaps. In Fig. 10 we show the overlaps 

As discussed in Section 6.2, the two 1~ 



for these states with operators {n, p} x D J=] 
hybrids are mixtures of spin-singlet and spin-triplet basis states and so they over lap with 
both the 7T x D^ =1 and p x D^ =1 operators. Therefore, for these states we plot tJ~Z% + Zp 
where Z n and Z p are the overlaps with these two different operators. The similar Z values 
suggest that these four states have a common structure, supporting their identification as 
members of a supermultiplet. 

The observed pattern of states in the supermultiplet, expected if a quark-antiquark 
pair in S-wave is coupled to a 1 H excited gluonic field, and its appearance at an energy 
scale ~ 1.2 GeV above the lightest conventional state are consistent with what was found 
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Figure 10. Overlaps, Z 1 of the D s states proposed to be members of the lightest hybrid supcrmul- 



tiplet with operators that are proportional to the field strength tensor, {tt, p} x D J=l 
into the relevant irreps. 



subduced 



in the light meson [37], light baryon [16] and charmonium [6] sectors, suggesting common 
physics. Those studies also discuss the first excited hybrid supermultiplet. We find candi- 
date positive-parity hybrid states in both the charm- light and charm-strange spectra ~ 1.5 
GeV above the ground state, the highest four positive-parity states in Figs. 7 and 8, but 
we do not robustly determine all the states in that energy region. 

We have given an interpretation of our results in terms of non-relativistic quark model 
multiplets. In the heavy-quark limit, where the charm quark is much heavier than the 
light and strange quarks, another classification scheme is convenient because the spin of 
the heavy quark decouples. The light degrees of freedom have total angular momentum 
j = L (g) s q , where s q = 1/2 is the spin of the light quark, and for L > this results in 
two possible values, j = \L ± 1/2|. Coupling with the spin of the heavy quark gives a 
pair of values for J for each value of j; this pair corresponds to a degenerate doublet of 
states in the heavy-quark limit, labelled by j p . For example, for S-wave (L = 0) states 
there is one doublet with j p = \ [(0, 1)~] and for P-wave [L = 1) states there are two 
doublets, | + [(0, 1) + ] and | + [(1,2) + ]. We discuss the heavy-quark limit in the context 
of mixing between spin-singlet and spin-triplet states in Section. 6.2. Although the charm 
quark is heavy and so these considerations may provide a useful guide, we do not expect 
it to be heavy enough for the heavy-quark limit to apply rigourously. This is observed in 
the spectrum where, for example, the 0~ and 1~ are not degenerate. 

Applying an analogous argument to hybrids, coupling the light quark degrees of free- 
dom with a 1 + excited gluonic field we have j = L(g)Sg(S> 1- For the lightest supermultiplet 
(with L = 0) this gives a j p = \ doublet with J p = (0, 1) _ and a j p = | doublet with 
J p = (1, 2) _ . In our spectra we note that, apart from the charm-light 2~ state, the (0, 1)~ 
and (1, 2)~ pairs of hybrids are more closely degenerate within a pair than between pairs, 
as suggested by the preceding argument. 

6.1 Comparison with experiment 

We now compare our results to the current experimental situation. Masses from the sum- 
mary tables of the PDG review [1] are shown in Figs. 7 and 8 where the pattern of the 
rather limited number of experimental states is seen to be in general qualitative agreement 
with our spectra. However, some comments on the quantitative differences are warranted. 
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As discussed in Section 5, states above threshold can have large hadronic widths and 
a conservative approach is to only consider our masses accurate up to this scale. An 
additional systematic uncertainty arises from having light quarks that are unphysically 
heavy leading to a pion mass of ~ 400 MeV (our strange quarks have approximately the 
correct mass). Below threshold we would expect this to be relatively unimportant in the 
charm-strange sector, but near or above thresholds involving mesons containing light quarks 
this could be significant. Another source of systematic uncertainty arises from having a 
finite lattice spacing and we discuss this in more detail below. We argued in Section 5 that, 
with a few exceptions, we see no significant finite-volume effects in our results. 

We obtain a value of 124(1) MeV for the charm-light S'-wave hyperfine splitting which 
is 18 MeV lower than the experimental value of 142.12±0.07 MeV for the neutral mesons [1]. 
We get 120(1) MeV for the analogous charm-strange splitting which deviates by 24 MeV 
from the experimental value of 143.8 ± 0.4 MeV. In our recent study of the charmonium 
spectrum [6] we found a short-fall of ~ 40 MeV in the S-wave hyperfine splitting. In- 
vestigating the dependence of the splittings on the spatial clover coefficient, c s , we found 
that an increase of this coefficient brought the hyperfine splittings in line with experi- 
ment. This suggested that the discrepancy was due to lattice discretisation effects and 
gave an approximate scale of ~ 40 MeV for the leading 0(a s ) systematic uncertainty. We 
performed an analogous test for the ground state S'-wave charm-light and charm-strange 
mesons and found a similar result; boosting the clover coefficient to c s = 2 from our tree- 
level tadpole-improved value (c s = 1.35) gave S-wave hyperfine splittings of 144(1) MeV 
in the charm-light sector and 139(1) MeV in the charm-strange sector, in line with the 
experimental values. This suggests that we can follow the arguments of Ref. [6] and assign 
a scale of ~ 20 MeV to our leading 0(a s ) systematic uncertainty. 

In the charm-light case, we find our P-wave states are heavier than the experimental 
equivalents; this could be due to the unphysically heavy mass of our light quarks and/or 
interaction with the nearby thresholds. In the charm-strange sector, two of our P-wave 
states are consistent with experiment but the other two states, expected to correspond to 
the enig matic L»* (2317) ± and L' s i(2460) ± , are significantly higher than their experimental 
counterparts. We note that the + and 1 + are very close to, respectively, the DK and D*K 
thresholds, and both the experimental and calculated + states lie the same distance from 
their appropriate thresholds. This may suggest that the unphysically heavy light quarks 
are a major contribution to the discrepancy. However, because of the interaction with the 
threshold, as discussed in Section 5, further study is required with multi-hadron operators 
included in our bases. 

We note that there have been experimental observations of additional higher-lying 
states; these are presented in the PDG's meson listings but omitted from the summary 
tables [1]. Because they are still awaiting experimental confirmation we do not include 
them in our plots. 

6.2 Mixing of spin-singlet and triplet states 

In QCD, the charm quark is significantly heavier than the strange and light (up, down) 
quarks and therefore St/(4) flavour symmetry is badly broken. In particular, charm-strange 
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J p 


~ (p — P2) ~ 7T ~ 7T2 


Heavy-quark limit 


c-1 


1 + 
2" 

1 _ (hybrid) 


6U.1(U.4J bz.6(U.2J o5.4(U.2J 
26.7(2.2) 22.2(3.7) 18.9(3.9) 
59.7(1.1) 68.4(0.8) 67.4(0.9) 


54. ( or 35.3 
50.8 or 39.2 


c-s 


1+ 
2 - 

1 _ (hybrid) 


60.9(0.6) 64.9(0.2) 66.4(0.4) 
64.9(1.9) 68.7(2.0) 70.9(1.8) 
59.9(1.7) 67.9(0.9) 67.3(0.9) 


54.7 or 35.3 

50.8 or 39.2 



Table 5. The absolute value of the mixing angles for the lightest pairs of 1 + , 2~ and hybrid 1~ 
states in the charm-light (c-1) and charm-strange (c-s) sectors. The angles extracted using different 
operators are presented; these are labelled by the gamma matrix structure with the derivative 
structures described in the text. Also shown are the mixing angles expected in the heavy-quark 
limit. The apparent difference between the charm-light and charm-strange 2 _ mixing angles is 
explained in the text. 



and charm-light mesons are not eigenstates of a generalisation of charge-conjugation. In 
contrast, SU(2) isospin symmetry is a reasonable approximation in the light quark sector 
and light mesons are eigenstates of G-parity to a good approximation 2 . The absence of 
such a symmetry for charm mesons means that S = 1 and S = quark-model states with 
J = L ( 3 Lj = l and Lj = l) can mix and probing this mixing can help us understand the 
internal dynamics and quantify flavour symmetry breaking. For flavourless mesons these 
would correspond to J PC = J ±=l: and J ±=F respectively. Using a two-state hypothesis and 
assuming energy-independent mixing (the states are close in energy), we can expand the 
pair of states A and B (we have chosen B to be the heavier state of our pair) in terms of 
spin-singlet and triplet basis states, 



\A) = +COS0 \ l L J=L )+ S \n9 \ 3 L J=L ) , 
\B) = -sme \ l Lj= L ) +cos0 \ 3 L J=L ) . 



3.1) 



In the non-relativistic limit, our operator (p — P2) x D 



states and 
(P 



{7r,7T 2 } x D 



J=L 



only overlaps onto 3 Lj = l 



erators, 
spin-trip 



p 2 ) x D 



[2] 



J=L 

and 



only overlap onto 1 -Lj=l states. There are analogous op- 



, which overlap onto respectively the 



{ir,ir 2 } x Df ] 
J=i I 1 Jj=i 

et and spin-singlet 1 _ hybrids. To determine the mixing angle 9, we take the 
ratio of the overlap factor of one of the aforementioned operators for state A to the overlap 
factor for state B; this ratio gives tan# or cot# depending on the operator used. 

Extracted mixing angles for the lightest pairs of P-wave (1 + ), D-w&ve (2~) and hybrid 
(1~) states are presented in Table 5 for each of the operators considered. Note that the 
overall sign (or in general the overall phase) of the overlaps for each state is not observable 
in our calculation and hence we determine the absolute value of the mixing angles. For 
each pair of states, the variation between mixing angles determined using the three different 



2 Our lattice calculations have degenerate up and down Wilson quarks so isospin is an exact symmetry. 
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operators gives an idea of the size of the systematic uncertainties. As well as the usual 
lattice systematics, these include uncertainties arising from this simple interpretation, a 
two-state hypothesis with energy-independent mixing, and from assuming that the energy 
difference between the states is sufficiently small such that the renormalisation factors for 
the operators do not vary substantially over this energy range. A relatively large variation 
is seen in the 1~ hybrid mixing angles; this would be expected because here the assumptions 
are less justified: these states lie higher in the spectrum, there are other states relatively 
close by with the same quantum numbers and there is a relatively large mass splitting 
within a pair. 

The P-wave ( 3 Pi - 1 P±) mixing angles are similar in the charm-strange and charm- 
light sectors, as are the mixing angles for the vector hybrids. For the D-waves ( 3 Z?2 - l D2), 
r _, rj 90° - 9 c —i because the dominantly 3 Z?2 state and the dominantly D2 state are 
almost degenerate, and the mass ordering of our states is reversed in the charm-strange 
sector compared to the charm- light sector. Although our strange quark mass is close 
to the physical value [8], our light quarks are unphysically heavy giving a pion mass of 
~ 400 MeV. Therefore SU{3) (strange- light) flavour symmetry is not badly broken (e.g. 
Mk/M-k is close to unity [12]) and we would expect to find similar mixing angles in the 
charm-light and charm-strange sectors. Because we have unphysically heavy light quarks, 
our determinations of mixing angles in the charm-strange sector might be expected to be 
closer to the physical values than those in the charm-light sector, but we note the other 
sources of systematic uncertainty discussed above. 

In Table 5 we also show the mixing angles expected in the heavy-quark limit where 
m c 3> m Uj d, s [3]. As discussed above, in this limit the heavy-quark spin decouples and 
states appear in doublets labelled by j p where j is the total angular momentum of the light 
degrees of freedom. For P-wave states there are two doublets \ [(0, 1)+] and § [(1,2)+]. 
The j p 1 + states (appropriate for the heavy-quark limit) are linear combinations of the 
singlet ( x Pi) and triplet ( 3 Pi) states (appropriate for the equal-mass limit) corresponding 
to a particular mixing angle. An analogous situation arises for the other pairs of states. 

Our mixing angles lie between zero mixing (the flavour symmetry limit, i.e. 0° or 
90°) and the heavy-quark limit values. This is consistent with what one expects because 
the charm quark mass is at an intermediate scale, larger than the light and strange quark 
masses but not large enough for the heavy-quark limit to apply. Quark model calculations 
and other determinations [2, 38-42] in general find significant mixing between the spin- 
singlet and spin-triplet states with, at least for the lightest 1 + pair, some deviation from 
the heavy quark limit. In comparison, because SU(3) flavour symmetry is not strongly 
broken our calculations, near zero mixing was found in the kaon sector [12]; we would 
expect the analogous mixing angles in the B meson system to be closer to the heavy-quark 
limit values. 

7 Summary 

We have computed extensive spectra of charm-light and charm-strange mesons using dis- 
tillation and the variational method combined with a large basis of carefully constructed 
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operators. Along with our spin identification scheme, these techniques have allowed us to 
extract a high number of states across all possible J p combinations, up to and including 
states of spin four. Our calculations were performed on two volumes, and apart from the 
exceptions mentioned earlier, we do not observe any significant volume dependence. 

Much of our spectra appears to follow the n 2S+1 Lj pattern expected by quark models 
but we also observe an excess of states that do not appear to fit into such a scheme. 
We have suggested that these states be interpreted as hybrid mesons due to their strong 
gluonic contribution. Their pattern ( J p = [(0, 1, 2)~, l - ]) and energy scale are consistent 
with that previously observed in both the light meson [37] and charmonia [6] sectors, and 
can be interpreted as a colour-octet quark-antiquark pair coupled to a 1 H chromomagnetic 
excitation. 

Using a non-relativistic interpretation of some of our operators, we extracted mixing 
angles between the spin-singlet and spin-triplet axial-vector 1 + states, 2~ D-wave states 
and 1~ hybrid states. We found mixing angles intermediate between zero mixing and 
the heavy quark limit values. Results in both the charm-light and charm-strange sectors 
are very similar, but we note that St/ (3) flavour symmetry breaking is suppressed on our 
ensembles as we have unphysically-heavy up and down quarks; Mk/M^ is close to unity. 
Nevertheless, we expect the flavour symmetry breaking in the charm-strange sector to be 
of the correct scale, bearing in mind the other systematic uncertainties discussed above. 

A full understanding of the charm sector requires us to study resonances and this is 
particularly relevant for many of the enigmatic charmed and charmonium-like states. We 
have calculated dispersion relations for the low-lying charm-light and charm-strange states, 
laying the foundations for a study of scattering states, along the lines of Refs. [26, 27], using 
carefully constructed multi-hadron operators and the Liischer methodology [28] with its 
extensions. 
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A Tables of results 



In Tables 6 and 7 we tabulate, respectively, the charm-light (D) and charm-strange (D s ) 
masses calculated on the 24 3 ensemble, as presented in Figs. 7 and 8. The scale has been 
set using the Q-baryon mass as described in Section 2; note that M n ~ 400 MeV in these 
computations. Half the calculated mass of the rj c has been subtracted in order to reduce 
the systematic uncertainty arising from any imprecision in tuning the bare charm-quark 
mass (see Section 2). 



J p 


M-M Vc /2 (MeV) 


Q- 


403(1) 


1155(7) 


1575(16) 


1811(14) 






1- 


527(1) 


1220(7) 


1361(5) 


1612(14) 


1704(19) 


1907(14) 


2 _ 


1377(5) 


1380(5) 


1856(15) 


1911(29) 


1944(19) 




3 _ 


1384(9) 


1990(17) 










4 _ 


1968(24) 


2028(23) 










0+ 


854(3) 


1505(11) 


1861(16) 








1+ 


959(3) 


992(3) 


1563(11) 


1565(11) 


1849(21) 


1919(13) 


2+ 


1024(3) 


1594(10) 


1707(5) 


1929(13) 






3+ 


1708(6) 


1718(8) 










4+ 


1720(14) 













Table 6. Summary of the D meson spectrum, as presented in Fig. 7, with statistical uncertainties 
shown. Note that M w ~ 400 MeV in these computations and more details are given in the text. 



J p 


M -M V J2 (MeV) 


Q- 


469(1) 


1186(8) 


1631(13) 1835(15) 






1- 


589(5) 


1267(6) 


1398(5) 1674(13) 


1771(16) 


1952(13) 


2~ 


1425(4) 


1429(4) 


1840(18) 1977(21) 


1986(15) 




3" 


1452(6) 


2030(14) 








4" 


2023(14) 


2027(16) 








0+ 


920(2) 


1543(10) 


1945(8) 






1+ 


1019(3) 


1055(5) 


1601(10) 1613(9) 


1895(16) 


1988(34) 


2+ 


1087(2) 


1651(10) 


1736(7) 1998(14) 






3+ 
4+ 


1749(7) 
1766(14) 


1767(7) 









Table 7. Summary of the D s meson spectrum, as presented in Fig. 8, with statistical uncertainties 
shown. Note that M„ ~ 400 MeV in these computations and more details are given in the text. 
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